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ABSTRACT 

Since the discovery of the first directly-imaged, planetary-mass object, 
2MASS 1207 b, several works have sought to explain a disparity between its 
observed temperature and luminosity. Given its known age, distance, and spec- 
tral type, 2MASS 1207 b is under-luminous by a factor of ~10 (~2.5 mags) 
when compared to standard models of brown-dwarf /giant-planet evolution. In 
this paper, we study three possible source s of 2MASS 1207 b's under-luminosity. 



First, we investigate iMohanty et al.l ( 120071 ) 's hypothesis that a near edge-on disk. 



comprising large, gray-extincting grains, might be responsible for 2MASS 1207 
b's under-luminosity. After radiative transfer modeling we conclude that the 
hypothesis is unlikely due to the lack of variability seen in multi-epoch photom- 
etry and unnecessary due to the increasing sample of under-luminous brown- 
dwarfs/giant-exoplanets that cannot be explained by an edge-on disk. Next, we 
test the analogous possibility that a spherical shell of dust, could explain 2MASS 
1207 b's under-luminosity. Models containing enough dust to create ~2.5 mags 
of extinction, placed at reasonable radii, are ruled out by our new Gemini/T- 
ReCS 8.7/im photometric upper-limit for 2MASS 1207 b. Finally, we investigate 
the possibility that 2MASS 1207 b is intrinsically cooler than the commonly used 
AMES-DUSTY fits to its spectr um, and thus it is no t , in f act, under- luminous. 



New, thick cloud model grids bv iMadhusudhan et al.l (120111 ) fit 2MASS 1207 b's 



l-10;^m SED well, but they do not quite fit its near-infrared spectrum. However, 
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we suggest that with some "tuning", they might be capable of simultaneously 
reproducing 2MASS 1207 b's spectral shape and luminosity. In this case, the 
whole class of young, under-luminous brown-dwarfs/giant-exoplanets might be 
explained by atmospheres that are able to suspend thick, dusty clouds in their 
photospheres at lower temperatures than field brown-dwarfs. 



1. Introduction 

Until recently, the vast majority of our knowledge regarding exoplanets has come from 
radial-velocity studies. Now with ground and space-based dedicated instruments, transit 
studies are discovering and characterizing the atmospheres of hot jupiters/neptunes at a 
rapid pace. The frontier for exoplanet studies is the direct detection and characterization of 
exoplanets, using high-contrast imagers/spectrographs. The emphasis in the direct detection 
field, thus far, has been on the discovery of planets. Now with a host of directly-detected 
planetary mass objects and candidates, we can proceed with multi-wavelength characteriza- 
tions of their atmospheres, the primary purpose of which, at least initially, is to determine 
the relationships between their masses, ages, radii, and spectra. 

In this paper we concentrate on one object, 2MA SS 1207 b, th e first directly de- 
tected, planetary-mass companion, which was found by IChauvin et al.l (|2004[ ) 0.78" from 
the brown dwarf 2MASSWJ 1207334-393254 (hereafter 2MASS 120 7 A). 2MASS 1207 A 
and its compan i on, 2M ASS 1207 b, were verified to be co-moving by IChauvin et al.l (120051 ) 
and ISong et al.l (120061). Since 2MASS 1207 A is known to be a member of the young TW 



Hya group (IGizisI 120021 ). IChauvin et all (I2Q0J assign 2MASS 1207 A and b an age of 8± 



Myr. Based on their measured JHKsL' fluxes, IChauvin et al.l ( 120041 ). estimate 2MASS 1207 



b's mass to be 5±2 Mjup with an effective temperature of 1250±200 K, assuming a distance 
of 70 pc. 



Several groups have refined the distance of 2MASS 1207 with the moving cluster method 



^The observations reported here were obtained at the Gemini Observatory, which is operated by the 
Association of Universities for Research in Astronomy, Inc., under a cooperative agreement with the NSF on 
behalf of the Gemini partnership: the National Science Foundation (United States), the Science and Tech- 
nology Facilities Council (United Kingdom), the National Research Council (Canada), CONICYT (Chile), 
the Australian Research Council (Australia), CNPq (Brazil), and CONICET (Argentina). 

^Based on observations collected at the European Southern Observatory, Chile (ESO Programmes 073. C- 
0469, 274.C-5036, 274.C-5057 and 083.C-0283). 

■^Based on observations collected with the SMARTS 1.3 meter, operated by the SMARTS consortium. 



- 3 - 



flBiUer fc Close 



2007: Gizis et al 



2007 



( Mamaiekl 12005: iMamaiek fc Meyerl l2007h and with trigonom etric parallax measurements 



Ducourant et al.l |2008| ) . For the remainder of this 



paper, we adopt the weighted average of the parallax measurements (52.8±1.0 pc) as the 
distance to 2MASS 1207. At this distance, the inferred luminosity of 2MASS 1207 b is 
considerably lower (~2.5 mags) than models predict given its age and spectral-type temper- 
ature]^ 



'he source of this under-luminosity is examined in great detail by iMohanty et al. 



(120071). With l o w-reso lution H/K spectra, and the use of accurate distance measurements. 



Mohantv et al.l (120071 ) f ind the temperatu re of 2MASS 1207 b to be 1600±100 K (using 
AMES-DUSTY models; lAUard et allboOll ) and the bolometric luminosity to be -4.72±0.14 
(log-f^). This luminosity is a factor of ~10 (~2.5 mags) less than the values predicted 



by AMES-DUSTY models. To explain their results, iMohanty et al.l ( 120071 ) invoke a "gray- 
extinction", edge-on disk model after arguing against distance/age errors, model errors, and 
interstellar extinction. A highly inclined disk around 2MASS 1207 b is, a priori, highly 
impr obable. Howeyer, 2MASS 1207 A is k nown to have a disk at a high (>60°) inclina- 
tion ( ISterzik et al.l 12004: IScholz et al.ll2005|). and disks aro und young binary stars tend to 



be close to aligned (jjensen et al. 



20041 : iMonin et al.l 120061 ) . Pe rhaps the best eyidenc e for 



the edge-on disk rnodel i s that other alternatives presented by IMohanty et al.l (120071 ) and 
Mamajek fc Meyerl (120071 ) are even more unlikely. 



The source of 2MASS 1207 b's apparent under-luminosity is important in the context 
of other low-mass objects, some of which (HD 203030 B, HN Peg B, and HR 8799 bcde) 
also appear to be less luminous than their temperatures sug gest, while others (GQ Lup b 



2007 



Mohanty et al. 



Marois et al. 



2007 



2008 



Close et al. 



2MASS 1207 A, AB Dor C and Beta Pic b) appear normal (iMetchev fc Hil 
Luhman et al.ll2007b 



Bowler et al.ll2010: Marois et al. 2010: McElwain et al 



enbrandl 



2006 



20071 . Bonnefoy et al. in prep). If 2MASS 1207 b's 



under-luminosity is indeed explained by an edge-on disk, other under-luminous low-mass 
objects might, at least individually, have their own reasons for being abnormal (such as 
inaccurate distance or age measurements). However, if 2MASS 1207 b's under-luminosity is 
not the result of an edge-on disk, it is very difficult to explain how an age/distance inaccuracy 
could explain its temperature/luminosity relationship, given its TWA cluster membership, 
established trigonometric parallax, and 2MASS 1207 A's normal temperature/luminosity re- 
lationship. Thus, alternative hypotheses should be considered that can explain 2MASS 1207 
b's under-luminosity in the context of other under-luminous brown-dwarfs/giant- planets. 



"^For the remainder o f this paper, we use the term " under-luminous" to describe this phenomenon. Other 
authors (see for example Metchev fc Hillenbrandl 2006 ) invert the problem using L — ^j^o-T^s, with a fixed 
(known) R, and describe Toff as "unusually cool" for its spectral type. 
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After presenting new observations and previously published data (Section [2]), we inves- 
tigate the source of 2MASS 1207 b's under-luminosity, in three parts, beginning with a test 
of the edge-on disk hypothesis (Section [3]). We first construct RADMC radiative transfer 
models of a hypothetical disk around 2MASS 1207 b. These models contribute to our dis- 
cussion of the expected variability of an edge-on disk, compared to observations of 2MASS 
1207 b. We also investigate the likelihood that a system like 2MASS 1207 b would be ori- 
ented on the sky nearly edge-on. Finally, we call attention to several other under-luminous 
brown-dwarf /giant-exoplanets that are difficult to explain with near edge-on disks. Analo- 
gous to the edge-on disk hypothesis, is the idea that an isotropic shell of dust could cause the 
same gray-extinction as an edge-on disk, without the geometric constraints (Section H]). We 
model this situation with the DUSTY radiative transfer code, and evaluate the models in 
the context of new, deep 8.7/im images from Gemini/T-ReCS, which provide a meaningful 
upper-limit on the mid-infrared brightness of 2MASS 1207 b for the first time. After study- 
ing the edge-on disk and dust shell hypotheses, both of which assume that 2MASS 1207 b is 
a Teff = 1600K brown-dwarf /giant-exoplanet with a low luminosity explained by extinction, 
we attempt to fit photometry and s pectroscopy of 2MASS 120 7 b with lower temperature 
thick-cloud model atmospheres from iMadhusudhan et al.l (120111 ) (Section E]). These models 
might be able to simultaneously explain 2MASS 1207 b's spectral shape and luminosity, in 
which case there would be no need to invoke either of the gray-extinct ing models. This solu- 
tion would also have the benefit of potentially explaining the whole class of under-luminous 
brown-dwarfs/giant exoplanets. 



2. Data and Observations 

We summarize the relevant published photometry/spectroscopy of 2MASS 1207 in Table 
[U We also present new Gemini/T-ReCS mid-infrared data, described in Section 12711 which 
allow us to put an upper-limit on the 8.7/im fiux of 2MASS 1207 b. We use these data 
to constrain disk models in Section [3l shell models in Section HI and brown- dwarf /giant- 
exoplanet model atmospheres in Section |3 

While photometry of 2MASS 1207 b has been published in a variety of near-infrared 
filters, multi-epoch data in any given filter rarely exists, the exception being NI CMOS pho 



tometry, which was published for two epochs in the F090M and F160W filters (ISong et al. 



20061 ). However, other archival data do exis t, that give more mu lti-epoch information than 



has previously been published. For example, IChauvin et al.l ( 120051 ) verify the common proper 



motion of 2MASS 1207 A and b, but do not publish the photometry from this data. We 
investigate 2MASS 1207 b's possible variability using archival VLT/NACO data in Section 
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We also explore 2MASS 1207 A's variability using SMA RTS /ANDICAM data that was 
previously published for parallax only (IBiller &: Close! 120071 ). This is relevant to much of the 
previously published photometry of 2MASS 1207 b, which was calculated relative to 2MASS 
1207 A and implicitly assumed 2MASS 1207 A was not variable. 



2.1. Gemini/T-ReCS 8.7/im Photometry 



We observed 2MASS 1207 in queue-mode on 2008 March 29 UT and in classical-mode 
on 2010 March 31 UT and 2010 April 1 UT at Gemini-South, with its highly sensitive 
Thermal- Region Camera Spectrograph (T-ReCS; iTelesco et al.l Il998l ). We used the Si-2 
filter {Xcentrai = 8.74/im), which is relatively insensitive to extinction/silicate absorption 
from a potential edge-on disk/shell, while being similarly se nsitive as the N-b and filteilfl, 
and mostly immune to variations in precipitable water vapor ( iMason et al.l 120081 ). Our data 
were taken in ~320 s on-source blocks, which corresponds to ~18 minutes clock-time when 
including chop-nod and other overheads. These long integration sequences are required to 
build up enough S/N to shift and add on 2MASS 1207 A. Combining the data in larger 
blocks appears to degrade our image quality due to guiding/nod-offset errors. 

We reduced our data with the custom T-ReCS IDL software MEFTOOLS v. 5.c|^, which 
allows the user to interactively display individual T-ReCS saveset frames and remove those 
with bad electronic artifacts/noise properties. We used this to discard approximately 2% of 
our frames and produce combined ~320 s chop-nod subtracted images. 

For each of these images , we fit 2MASS 1 207 A with a 2D Gaussian ellipsoid using the 
IDL software suite MPFIT (IMarkwardtl 120091 ). We discard images where the fit centroid 
error on 2MASS 1207 A is >0.5 pixels (0.045"). This keeps our final combined image free of 
degradation from shift and add errors (T-ReCS' image quality is typically ~0.3" FWHM). 
Discarding frames with large centroid errors also has the benefit of selecting the frames with 
the best image quality (FWHM image quality cannot be accurately measured because of the 
low S/N on 2MASS 1207 A in a single frame). Our selection leaves us with 15 ~320 s blocks, 
which are then weighted by the S/N of the Gaussian ellipsoid fit, and combined. While 14 
out of 25 blocks from 2010 March 31 UT were usable, only one out of the 6 blocks from 
2008 March 29 UT and none of the blocks from 2010 April 1 UT passed our centroid error 
cut. The image quality from 2010 April 1 UT appears to have been significantly degraded 



http : //www . gemini . edu/sciops/?q=node/10085 



http : //www . j im-debuizer . net/ research/ 
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by a nod-return/guiding error. A summary of our observations and weather conditions is 
presented in Table [21 Our final combined image is 4749 s on-source (9498 s open-shutter, 
including chop-nod subtraction), and is shown in Figured! 

Figure [1] shows 2MASS 1207 A at the center of the two red circles, and a green circle 



at th e near-IR determined position of 2MASS 1207 b (sep=0.773", PA=125.37°; ISong et al. 



2006|). The measured FWHM of 2MASS 1207 A is 0.30", so the core of the 2MASS 1207 
A point-spread-function (PSF) should have a negligible affect on the measured flux at the 
position of 2MASS 1207 b. However, there does appear to be a slight increase in background 
at a separations <1", which is probably the result of 2MASS 1207 A's PSF seeing halo. We 
subtract off this small halo contribution by doing a median average of all pixels within a 2.5 
pixel radius of the separation of 2MASS 1207 b, and at all position angles except for the 
aperture that we use for 2MASS 1207 b. 

We used data from 2010 Mar 31 UT to p erform an absolute flux calibration on 2MASS 
1207 A using the mid-IR standard HR4888 JCohen et al.lll999[ l. which is 16.011 Jy in the 
Si-2 filtei0. The best 4-hour period of 2MASS 1207 data (which included 9 of our 15 usable 
frames) was averaged, and compared to images of HR4888 taken immediately after (HR4888 
was found to vary by ~5% over the 3 calibration images). We use a 15 pixel (1.35") aperture 
and find that 2MASS 1207 is 0.000343 times as bright as HR4888 in the Si-2 filter, which 
corresponds to an absolute flux of 5.49 mJy. We assign an 8 % error to this measurement 
(0.44 mJy), which is conservative given the fidelity of the Si-2 filter and the night's good 
photometr ic quality. This valu e is consistent with published values of 5.6 mJy ±1 mJy in the 



Si-2 fi lter (jSterzik et al.ll2004j). 5.74 mJy in the slightly bluer Spitzer 8/im filter ( iRiaz et al. 
2006), Spitzer spectroscopy (jMorrow et aL 2008), and roughly consistent with our own mea 



surements using the 2008 Mar 29 UT data (4.3 mJy ±0.5 mJy), which suggests that 2MASS 
1207 A is probably not wildly variable at these wavelengths, despite the presence of a near 
edge-on disk. 

Using 2MASS 1207 A as a PSF, we do a best-fit to determine how muc h fiux is at 
the near-IR determined position of 2MASS 1207 b (sep=0.773", PA=125.37°; ISong et al. 



20061 ). We assume Gaussian error bars (equal on all pixels for the background limited case), 
and fit within a 5 pixel diameter aperture (see Figure [T]) to avoid contamination by 2MASS 
1207 A. We do a Monte Carlo evaluation of the background fiuctuations by doing a similar 
aperture measurements at 10000 randomly chosen position angles and separations > 1.5" 
from 2MASS 1207 A to avoid residual halo contributions. For each Monte Carlo trial, we 



'http : //www . gemini . edu/ sciops/instruments/midir-resources/imaging-calibrations/michelle-std-f luxes 

which is calculated for Gemini-N/Michelle, but the Si-2 filters of Michelle and T-ReCs are similar 
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subtract the flux of a randomly chosen background aperture from the measured flux at the 
position of 2MASS 1207 b. This gives us an estimate, for each trial, of what the true flux of 
2MASS 1207 b might be. We discard all negative flux measurements based on the Bayesian 
prior that 2MASS 1207 b's flux must be positive, and use the remaining cases to construct a 
cumulative distribution function, which calculates the probability that 2MASS 1207 b's flux 
is less than a given value (but greater than 0). We include an 8% absolute calibration error 
(see the previous paragraph), although increasing this error up to ~20% has a negligible 
effect on our results. We flnd with 50% confidence (i.e. median) that 2MASS 1207 b's flux 
is < 0.26 mJy, 84% confldence (la) that 2MASS 1207 b's flux is < 0.48 mJy and ~99.9% 
confidence (3a) that 2MASS 1207 b's flux is < 0.92 mJy. 



2.2. VLT/NACO Near-Infrared Photometry 

2MASS 1207 A and b have been spatially resolved by VLT/NACO at several epochs. 
Some of these data have been published as astrometry only, and some have not been pub- 
lished at all. Here we use archival data to investigate the possibility that 2MASS 1207 b is 
photometrically variable. 

2MASS 1207 was observed 6 times in the Ks filter between 2004 April 27 UT and 2009 
June 30 UT. Shorter wavelength filters were also used at multiple epochs, but generally on 
the same nights as Ks, so variability information from these filters is redundant. Shorter 
wavelength filters also suffer from lower signal-to-noise and a worse adaptive optics correction. 
Therefore we only focus on the Ks-band data here. 

Of the 6 epochs of Ks data, 2 dates (2005 Mar. 30 and 2005 May 18) suffered from 
poor seeing so that we could not produce reliable resolved photometry of 2MAS S 1207. The 



other four epochs were reduced using the custom IRAF pipeline described in IClose et al. 



(120031 ). We performed PSF- fitting photometry on 2MAS S 1207 A and b, and a background 



star (2MASS 12073400-3932586), using IRAF /daophot fIStetson] 119871 ) and allstar. Seeing 
conditions on 2009 June 30 were marginal, so the data were reduced by aperture photom- 
etry and assigned conservative errors. A summary of our observations and photometry is 
presented in Table [31 Ignoring the low-quality 2009 June 30 data, we find that compared to 
the background star, 2MASS 1207 A is variable by, at most, a =0.05 mag (0.09±0.03 mag 
peak-to-peak), which is similar to our optical SMARTS/ ANDICAM results in Section 1^751 
2MASS 1207 b is variable by, at most, a =0.12 mag (0.23±0.07 mag peak-to-peak), which 
is not abnormally high for a young object with star-spot and/or accretion variability. 
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2.3. SMARTS/ANDICAM I-band Photometry 



The combined system, 2MASS 1207 A and b was observed by iBiUer fc Close! (120071 ) 
between January 2006 and April 2007 with SMARTS/ANDICAM to measure the parallax 
distance to 2MASS 1207. Here we present the photometric results of this data set, which 
provide 53 nights of I-band photo metry, taken raonthly , in groups of several nights. Ob- 
servational details are described in lBiller fc Close! (12007! ). Special care was taken to always 
position the target on the same CCD pixel, and all data were taken near transit to aid the 
interpretation of the parallax data. Relative photometry between 2MASS 1207 and the 12 
brightest stars in the field gives an average variability of cr = 0.03 mag. Given that the 
differe ntial magnitude between 2MASS 1207 A and b is 7.8 in the similar NICMOS F090M 
filter (ISong et al.ll2006l ). this variability must be completely dominated by 2MASS 1207 A's 
intrinsic variability. The variability is small er than the error bar s of published relative pho- 
tometry of 2MASS 1207 b (see for example IChauvin et al.l 12004 ). which legitimizes 2MASS 
1207 A's use as a photometric calibrator for 2MASS 1207 b (although as we describe in 
Section 12. 2^ the nearby background star 2MASS 12073400-3932586 is a better calibrator for 
variability studies). 

Variations in the I-band magnitude of 2MASS 1207 with respect to its average I-band 
value are plotted in Figure |2J Error bars are calculated empirically from the multiple frames 
taken on each individual night. 



3. The Near Edge-on Disk Hypothesis 



A possible explanation for 2MASS 1207 b's under-luminosity is th at it is being partially 
extincted by a near e dge-on disk of large (gray-extincting) dust grains (IMohanty et al.!l2007l : 
Patience et al.!l2010l ). Such a disk must produce ~2.5 mags of extinction at J-band, with 
nearly gray extinction. The disk must create only moderate photometric variability. The 
disk may not have an extremely unlikely viewing angle/geometry. And the disk properties 
should be consistent with our current picture of disk evolution. Additionally, any dust around 
2MASS 1207 b will have its emission constrained by our 8.7 fim photometry. In this section, 
we discuss the viability of the near edge-on disk hypothesis, with respect to these conditions. 



3.1. RADMC Disk Models 



We used the RADMC radiative transfer code (iDuUemond fc DominikI 120041 ) and RAY- 
TRACE, a post-processing tool, to calculate the emerging spectral energy distribution (SED) 
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of a hypothetical disk around 2MASS 1207 b. The primary input parameters of the RADMC/RAYTRACE 
codes are the SED and mass (Mgtar) of a central source, and the inner radius (-Rm), outer 
radius (-Rdisk)? inclination (z), and mass (M^isk) of its disk. Other relevant parameters in- 
clude the disk's vertical flaring geometry (described below), and a dust grain distribution 
with associated optical properties. 

As we will show later in this section, our Gemini/T-ReCS 8.7/im photometry upper-limit 
is not sensitive enough to detect the presence of a disk around 2MASS 1207 b. Consequently, 
there is no data available at wavelengths longer than L-band to distinguish between different 
possible disk models. However, 2MASS 1207 A has a considerable amount of data at longer 
wavelengths. As such, our strategy is to model the disk around 2MASS 1207 A first, and 
scale the model for use around 2MASS 1207 b. 

For 2MASS 1207 A, we adopted a T,ff=2600 K, Ra=0 25 Rq, log{g)=A and Ma=0.025 
Mq AMES-DUSTY synthetic sp ectrum (lAUard et al.ll2001l) for the central source, similar to 
the best-fit model determined by lMohanty et al. fl2007h ~ For t he disk parameters, we assume 
-Rin=3.3 -Ra (the dust sublimation radius for 2MASS 1207 A; 



Morrow et al.ll2008h . i?disk=20 



AU (~0.46xprojected binary separation: lArtymowicz fc Lubowlll994j ). and Mdisk=0.01 Ma. 
We note the precise values of -Rdisk and Mdisk have a negligible effect on the SED for A <20/im. 

We vary the model's dust grain distribution, vertical geometry, and disk inclination to 
determine a best-fit to the data desc ribed in Section El For the dust grain distribution we 
applied a single olivine dust species (jJager et al.l 120031 ) and varied the minimum diameter 
(omin) of the grain s between 0.1 /im a nd lO/im. We assumed a dust grain size power-law 



with n{a) oc a~^'^ ( iMathis et al.l 119771 ). and a maximum dust grain size of 1 mm. For the 
vertical geometry, we ran two sets of models: a fully- flared model (hereafter, "flared"), where 
RADMC calculates the disk-height self-consistenty (assuming hydrostatic equilibrium and 
dust-gas coupling), and a set of flatter (hereafter, "flat") models where the disk scale heights 
are parameterized. The flat disk models are assumed to follow a power-law description with 
^ oc r^/'^ and outer disk pressure scale heights (Hp/r) ranging from 0.04 to 0.10 at -Rdisk- 
The final free parameter, inclination, is varied between 5° and 90° to determine a best-fit. 

The flared disk models do not fit the observed SED of 2MASS 1207 A. While the models 
can fit the near-IR observations (assuming inclinations <65°), they predict significantly 
higher fluxes than are observed in the mid-infrared (See Figure [3]). The best-fit flat disk 



model {Hp/r=OM, 



-7fim and z=71°), also shown in Figure [3], provides an adequate 



fit to all of the 2MASS 1207 A data, although other parameter combinations also result 
in acceptable fits. Exploring the degeneracies in these parameters, we find that the disk 
around 2MASS 1207 A has 0.06< H^/r <0.08, a^i„ >5/xm, a r id 70° < i < 75°. The best-fit 
flat models are consistent with the results of iMorrow et al.l (120081 ). Both models suggest 
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advanced dust processing and dust settling in the disk around 2MASS 1207 A. 



To model the hypothetical disk around 2MASS 1207 b, we scaled down the flared and flat 
disk models of 2MASS 1207 A. For the central source we use an AMES-DUSTY model with 



Teff=l,600 K, i?A=0.16 Rq and log{g)=4:.5 (jPatience et al.ll2010l ). We used the relevant disk 
parameters from 2MASS 1207 A: Hp/r=OM, Mdisk=0.01 Mb, Rn=1 14: (the dust subli- 
matio n radius) and -Rdisk=10 AU (~0.2xprojected binary separation: lArtymowicz fc Lubow 
19941 ). Our best-fit fiared and fiat models are shown for 4 different values of Omin in Figure 



m The lack of (A >10/im) mid- infrared data makes it impossible to distinguish between the 
fiared and fiat disk models, although we note that disks aroun d low-mass stars/b rown dwarfs 
are found to be fiatter than disks around higher- mass stars (ISziics et al.ll2010l ). We varied 
ajnin from O.l/xm to 5/im and find that the amin=0.1/im fiat-disk models cannot simultane- 
ously fit the Kg and L', and NICMOS observations. Models with a^i^ >0.5yum fit the data 
adequately, although the disk models with amin = l/wm or 5/im do not fit the data quite as 
well as the Omin = 0.5/im due to the redness observed between L' and the shorter wavelength 
data. In each case, the models have been fit by varying inclination, which is constrained by 
our requirement that the disk extinct the central source by ~2.5 mags at J-band in order to 
fit 2MASS 1207 b's near-infrared photometry. Figure O shows the amin =0.5/im fiat model 
from Figure m but with varying inclination. This model is consistent with the measured data 
for an inclination range of 80.6° < i < 80.9°. Different disk models have similarly tightly 
constrained inclinations, although the actual ranges vary between models. For a given model, 
the inclination of the disk (or non-axisymmetric structure) cannot significantly vary in time 
without producing strong variability (see Section | 



For the remainder of this paper we use the amin =0.5/im models for 2MASS 1207 b. We 
continue to discuss the fiared and fiat models, although we prefer the fiat models due to the 
presence of one around 2MASS 1207 A, an d more generally, because of the ubiquity of fiat 
disks around low-mass stars/brown-dwarfs (jSziics et al.ll2010l ). 



3.2. Variability in Near Edge-on Disks 

The lack of near-infrared variability exhibited by 2MASS 1207 b would be unusual for a 
system with a partially extincting, near edge-on disk. Young systems with edge-on or near 
edge-on disks have been shown to exhibit variability over a variety of masses (brown-dwarfs, 
T Tauri stars, Herbig Ae-Be stars). This variability is generally the result of eclipsing non- 
axisymmetric structures in Keplerian orbit, such as warps, hydrostatic fiuctuations, gaps, 
spirals, and dumpiness. The best studied edge-on variables are UX Ori type stars, which dis- 
play frequent eclipse-like extinction events caused by hydrodynamic fiuctuations along their 
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puffed-up inner-rims. However, the UX Ori phenomenon is thought to o nly exist in com- 



plete ly self-shadowed disks, which are usually around Herbig Ae-Be stars (jPuUemond et al. 



20031 ). In lower-mass T-Tauri stars and brown dwa rfs, disk flaring causes the outer-parts of 
the disk to be the dominant source of extinction (jPuUemond et al.l l2003l ) . The structural 
sources of variability in the outer parts of disks are less well understood than inner-rim vari- 
ability. A far-away companion, mis-aligned wit h the plane of the disk (i.e. 2MASS 1207 A) 
can cause warps in the outer-parts of the disk (IFragner fc Nelsonll2010l ). but in this case the 
warp will orbit with the period of the binary, which is much longer than the orbital period of 
the disk, and is unlikely to cause short-term variability. Other source s of non-axisymmetric 
struc t ure have been stud ied with regard to specific perturbations (IFlaherty &: MuzeroUe 
2010l : iFlaherty et al.l 1201 ll ). but a general theory does not yet exist. 



Empirically, we know that edge-on and near edge-on disks that extinct their host star 
exhibit variability over many timescales. Famous examples include t he T Tauri stars AA 



Tau, which exhibits ~1 mag variability on periods of days/weeks (iBouvier et al.l 11999 
20031 ) , HH 30, which shows dramatic variability in spa tially resolved scattered light im 



ages (IStapelfeldt et al.lll999l : IWatson fc Stapelfeldtl 120071). T Tau Sa, which varies in bright 



ness by several magnitudes on timescales of years (ISkemer et al.l l2008t Ivan Boekel et al. 
2010l). and UY Aur B, which became fainter by several magnitudes over a ~50 year period 



(IJov &: van Biesbroeckl Il944j : iHerbst et al.lll995l ). A systematic analysis of NGC 2264 by 
Alencar et al.l (|2010[ ) revealed that 28%±6% of classical T Tauri s tars exhibit AA Tau- like 
variability with amplitudes as large as 137% over a 23 day period. lAlencar et al.l (120 lOf ) do 
not publish the variability amplitude of all of their AA Tau-like systems, so we do not know 
the lower-limit amplitude. However, their result systematically confirms that AA Tau-like 
variability is common, and probably synonymous with near edge-on/edge-on disks. 

The physical mechanisms leading to disk asymmetries vary as a function of spectral 
type and accretion rate and thus the well-studied disks described in the previous para- 
graph may not be ana l ogous to brown dwarfs. Although exara ples are sparse, several groups 
( iLuhman et al.ll2007at IScholz et al.ll2009l : iLooper et al.ll2010bl Jal) report evidence for variable 
occultation of brown dwarfs from edge-on or near edge-on disks. These cases provide em- 
pirical examples, but unbiased, systematic variability studies of near edge-on/edge-on disks 
around brown dwarfs are rare. 

From a modeling perspective, a small warp of just 1 degree (which we approximate as 
changing the inclination of our RADMC disk by 1 degree) would cause variability of am- 
plitude 2.20 mags in the NICMOS F160W filter (as well as 2.13 in the NICMOS F090W 
filter, and 2.16 in Ks-band). This is inconsistent with measured variability amplitudes of 
0.03±0.03 in F160W, 0.24±0.5 mag in F090W, and 0.23±0.07 in Ks-band as well as the 
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general stability of the multi-epoch, multi-wavelength photometry and spectroscopy sum- 
marized in Table [TJ With respect to the variable browii dwa rfs discussed above, w hich we 
caution is not an unbiased sample, iLooper et al.l (l2010bl ) and iLooper et al.l (l2010al ) find an 
amplitude of variabi lity of 1.3 mags and 0.9 mags at Ks for TWA 30A and TWA SOB over 
a period of months, IScholz et al.l (j2009f ) find an amplitude of variability of 0.45 rn ags at K 
for their "Object #2" over a period of about a week and iLuhman et al.l (l2007al ) find an 
amplitude of variability of ~0.5 mags at Ks for 2MASS 0438-1-2611 over 2 epochs separated 
by 6 years. All of these objects are substantially more variable than 2MASS 1207 b. 

As a result, we find the existence of a near edge-on disk around 2MASS 1207 b to be 
incompatible with 2MASS 1207 b's measured photometric stability. However, we caution 
that this conclusion is based on an incomplete knowledge of the structure and variability 
of extincting disks around brown dwarfs. Further studies that catalog lightcurves of many 
young low-mass stars (such as the ongoing Spitzer/YSOVAR program, and future surveys, 
such as PAN-STARRS and LSST) will better constrain the expected variability of extincting 
disks around young brown dwarfs. 



3.3. The Expected Frequency of Partially Extincting Disks Around Young 

Brown Dwarfs 

A priori, it would be unlikely (and perhaps unfortunate) that 2MASS 1207 b, a system 
which is currently unique as the closest, youngest, extremely low-mass object, is inclined to 
our line of sight in a near edge-on configuration where it is "partially" extincted. In this 
section, we investigate the probability that a brown dwarf system is in a configuration such 
that its disk "partially" extincts the central source. We use our RADMC disk models for 
2MASS 1207 b and calculate the extinction produced by the hypothetical disk when viewed 
at different inclinations varying from to 90 degrees. 

In Figure in] we show the H-band extinction predicted by the models as a function of disk 
inclination, for both fiat and fiared geometries. We note that there are 3 distinct regions 
of each plot: an unextincted region, a partially extincted region where the extinction rises 
sharply, and a fully extincted region where the extinction levels off. Physically, the partially 
extincted region begins when the outer edge of the disk moves into our line of sight towards 
the brown dwarf. As inclination continues to increase, the extinction towards the brown 
dwarf rises until the fully extincted region, where scattered light (which is assumed to be 
isotropic in RADMC) caps extinction even as opacity continues to rise. 

Based on this analysis, and assuming an isotropic distribution of disk inclinations, there 



- 13 - 



is a 81.5% probability that a brown dwarf with our flat disk model is unextincted [Ah =0-0.5 
mags), a 6.4% probability that it is partially extincted {Ah =0.5-6.0 mags), and a 12.1% 
probability that it is fully extincted [Ah =6.0-7.1 mags). For the flared disk model, there 
is a 51.0% probability that it is unextincted {A^ =0-0.5 mags), a 14.6% probability that it 
is partially extincted {Ah =0.5-3 mags), and a 34.4% probability that it is fully extincted 
{Ah =3-4 mags). For brown dwarfs, the fiat disk model is more accurate than the flared 
disk model (see Section |3TT]) . Although the exact analysis depends strongly on disk flaring, 
the two cases presented here bracket the likely geometries. 

Assuming 2MASS 1207 b has a disk, but ignoring its under-luminosity as prior evidence 
for why it should have a near edge-on disk, the probability that 2MASS 1207 b's disk is 
in the inclination range that causes partial extinction is ~6.4%. This also neglects the fact 



(Jensen et al. 


2004; 


Monin et al. 


2006) 



object, and because we do not yet understand the likelihood that alternative effects are 
causing its under-luminosity, we cannot use the fact that 2MASS 1207 b's hypothetical disk 
is at an unlikely inclination to reject the disk's existenc^. However, when more objects like 
2MASS 1207 b (i.e. young, low-mass brown dwarf s/exoplanets) are found by surveys like 
WISE, our models predict that 81.5% will be unextincted, ~6.4% will be under-luminous 
by 0.5-6.0 mags, and 12.1% will be under-luminous by 6.0-7.1 mags. A higher percentage in 
the middle category would imply that disk inclination cannot explain the under-luminosity 
phenomenon. 



3.4. Other Under-luminous Brown-Dwarfs/Exoplanets 

While the previous section predicts the percentage of objects in 2MASS 1207 b's class 
that should appear under-luminous (partially extincted) due to the presence of a near edge- 
on disk, there are no other objects quite like 2MASS 1207 b in terms of its youth and low 
mass. However, there are other known under-luminous brown dwarfs and giant planets (e.g. 
HD 203030 B, HN Peg B and HR 8799 bcde), which are known to be older than 2MASS 
1207 b. Since the gas-rich, geometrically thick disks around these older systems will have 



R-w-\ n • T-» 7 ii 7-)/ 1 I J 1 ■ \ Plunder — luminous\edqe—Qn)P(edqe— on) ^T^ 

i'ollowiiia; Baves tneorem, Fyedqe — onmnaer ~ Luminous) — — 57 — — ', '-^ — -. io 

c> J 1 \ if I J Plunder— luminous). 

calculate the probability that a brown dwarf has a near edge-on disk, given that it is under-luminous, we 
need to know the probability that brown dwarfs are under- luminous (the denominator term). Because we 
do not understand the likelihood that alternative effects are causing the under-luminosity, and because we 
do not know what percentage of brown dwarfs are under-luminous, we do not know P[under — luminous) , 
so we cannot evaluate P{edge — on\under — luminous) 



- 14 - 



dissip ated (and in the case of HR 8799, the disks would more hkely be face-on; iMarois et al. 
20081 ). the near edge-on disk hypothesis proposed for 2MASS 1207 b would not be applicable. 
In the previous sections, we have stated why we find the 2MASS 1207 b near edge-on disk 
hypothesis to be unlikely. A further reason is that there are other under-luminous brown 
dwarfs that cannot be explained by a near edge-on disk, and thus, there is no reason to 
invoke one to explain the behavior of 2MASS 1207 b. 

Are there other patterns that might be shared between these under-luminous brown 
dwarfs? iMetchev et al.l (120091 ) note that all of the known und er-luminous objects are young 
L/T transition brown-dwarfs (although iPatience et al.l (|2010[ ) describe 2MASS 1207 b as an 
early L-type), and this pattern (sometimes parsed as the gravity dependen c e of the L/T 
transi ti on) has been studied b y numerous groups , inclu ding iLiu et al.l (2008). IDupuv et al. 



J2009h . lstephens et al.l ( l2009[ ) . Isaumon fc MarlevI toQ^ and lSowler et al.l J2OI0I). As brown 



dwar fs cool from L-type to T-type, their dusty clouds are thought to dissipate (IBurgasser et al. 
2OO2I ) and settle below the brown dwarf photosphere. This pattern implies that non- 
equilibrium chemistry and dust-cloud physics might play a crucial role in explaining these 
systems' apparent under-luminosities. 

A possib le explanatio n for th e un der-luminosity of the HR 879 9 planets has recently been 
proposed by Currie et al. ( 2011 ) and Madhusudhan et al. (j2011 ). who demonstrate that a 
thick cloud atmosphere can reproduce the observed photometry of the HR 8799 planets. The 
thick cloud models might be applicable to 2MASS 1207 b as well, and we investigate this 
possibility in Section 



4. The Dust Shell Hypothesis 

The problems in invoking a near edge-on disk to explain the under-luminosity of 2MASS 
1207 b stem from the fact that near edge-on and edge-on disks tend to produce more variabil- 
ity than is detected for 2MASS 1207 b, along with the fact that such a configuration is a priori 
unlikely and cannot be used to explain other under-luminous brown dwarfs. In this section, 
we model the system with a shell of dust, which should provide the extinction of a near 
edge-on disk, without the geometric effects that we have previously d eemed unlikely. Such a 



shell may have been discovered around the brown dwarf, G 196-3 B ( IZapatero Osorio et al. 



2OIOI). 
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We use the dust shell radiative transfer code, DUSTYIjI^ (llvezic et al.lll999l ). to model 
dust shells at a variety of distances from the central brown dwarf. In all cases, the central 
source is the AMES-DUSTY^° brown dwarf model atmosphere (Tefj=l,600 K, i?b=0.16 Rq 
and log{g)=4:.5) used as a model of 2MASS 1207 b in Section IXTl We also use the same dust 
grain size distribution and optical properties as described in I3.H with a minimum grain-size 
of amin=l/im (although varying this does not significantly change our results). For each 
shell, we fix the J-band extinction to be 2.5 magnitudes (the observed under-luminosity of 
2MASS 1207 b). The shells are geometrically thin, and placed at i?b, 5-Rb, 10-Rb, 20i?b, 
40i?b, and 80-Rb- Since DUSTY is a 1-D code, we are forced to choose a slab geometry 
(which we use for the i?b niodel), and a shell geometry (which we use for all others), where 
the central source is assumed to be a point-source. The difference between using a slab 
and point-like approximation is small compared to the bulk SED shape, so we consider the 
point-like approximation reasonable for our purposes. Another limitation of this modeling 
approach is that a shell of warm dust right near the surface of the brown dwarf would heat 
the brown dwarf and change the output spectrum. Thus our closest-in shells, while useful for 
modeling the mid- infrared output, are somewhat non-physical, especially in the near-infared, 
where self-consistent models would be necessary to produce the correct spectral features. 

The results of our shell models are presented in Figure [3 Dust shells inside of 20i?b are 
ruled out by our Gemini/T-ReCS 8.7/im +3a upper-limit, while homogenous shells outside 
of this radius are not observed for objects of 2MASS 1207 b's age and mass. 



5. The Thick Cloud Atmosphere Hypothesis 



JHK spectroscopy of 2MASS 1207 b shows its features to be consistent with a 1600 
K AMES- DUSTY model atmosphere, but its lumino sity is more consistent with a ~1000 
K model (IMohanty et al.l 120071 : iPatience et al.l l2010l ). Typically, dust is thought to set - 
tle/condense below the photosphere for atmospheres < 1200-1400 K (jSaumon fc Marleyll2008l ). 
so ~1000 K brown-dwarfs/giant- planets are not expected to exhibit significant dust opacity. 
However, models with cloud structures that are suspended above the photosphere in cooler 



^User Manual for DUSTY, University of Kentucky Internal Report, accessible at 
|http : // www . pa . uky . edu/~moshe/ dusty| 

^^Regrettably, the DUSTY dust-shell radiative transfer code shares a name with the AMES-DUSTY 
brown-dwarf model atmospheres we use in this paper. So to be clear, in this section, we are enshrouding an 
AMES-DUSTY brown-dwarf model atmosphere with a DUSTY dust-shell. Throughout this paper, we refer 
to the brown-dwarf models atmosphere's "AMES-DUSTY" (although they are often referred to by other 
authors as DUSTY) , and the dust-shell modeling software as "DUSTY" . 
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atmospheres might be able to explain the spectral features of 2MASS 1207 b at the low 
temperatures indicated by its luminosity. 



Recently, ICurrie et al.l ( 120111 ) and lMadhusudhan et al.l (120111 ) suggested that thick cloud 
atmospheres could explain the red colors and under- luminosity of the HR 8799 planets. While 
their models do an adequate job reproducing the planets' photometry, a lack of detailed 
spectroscopic data diminishes the significance of their test. In this section, we attempt to 
fit the available 2MASS 1 207 b data (see Table [T]) with the thick cloud atmospheres of 
Madhusudhan et al.l (l201l[ ). Medium-resolution, near-infared spectroscopy of 2MASS 1207 



b JPatience et al.lboioh provides the strongest test yet of these models. 



The iMadhusudhan et al.l (120111 ) models (see also iBurrows et al.l 120061 ) parameterize 
cloud thickness from thickest ( "Model-A" , where the clouds extend to the top of the atmo- 
sphere), to thin ("Model-E") with two intermediate specifications ( "Model- AE" and "Model- 
AEE"). Model gridfl also vary Tea, ^og{g), metallicity, cloud composition and modal grain 
size. The A-models maintain a constant mixing ratio with decreasing pressure, and obser- 
vationally are chara cterized by red near-i n frared colors. Similar colo rs are p roduced by the 
fsed = 1 models of lAckerman fc Marleyl (120011 ). ISaumon fc Marleyl (120081 ) and references 
therein, which might provide a similar explanation for 2MASS 1207 b's appearance. 



We plot a representative sample of the IMadhusudhan et al.l (120111 ) models, as well as a 
best-fit, scaled AMES-DUSTY model (jPatience et al.ll2010l ) against the observed photometry 
(Figure IHD and spectroscopy (Figure [H]) of 2MASS 1207 b. In both figures, frame (a) shows 
the best-fit AMES-DUSTY model, frames (b), (c) and (d) show a sample of A-models (thick 
cloud), and frames (e) and (f) show a sampl e of A E-models (intermediate between thick and 
thin clouds). All of the IMadhusudhan et al.l (120111 ) models we show assume log(g)=4.0, solar 
metalhcity, and a forsterite cloud composition. 



As described in IPatience et all tom . the 1600 K, log(g)=4.5 AMES-DUSTY fits the 
spectral shape of 2MASS 1207 b adequately. However, it can only fit the luminosity of 
2MASS 1207 b by scaling the model radius to 0.052 Rq, which is significantly smaller than 
the ~0.16 -RfT) radius predicted by evolution ary mod els for an object of 2MASS 1207 b's age 
JSurrows et all Il997t Ichabrier et aD I2OO9I ). Both IMadhusudhan et al.l J201ll )'s A-models 
(thick cloud) and AE-models (intermediate) are able to approximately fit the overall SED 
of 2MASS 1207 b without an unphy sical radius scaling ( the models simply assume a radius 
based on the evolutionary models of iBurrows et al.lll997l ). 



The JHK spectroscopy of 2MASS 1207 b provides a stronger test of the IMadhusudhan et al 



available at http : //www . astro . princeton . edu/ -burrows /8799/8799 . html 
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( 120111 ) models than the SED photometry. Of the AE models, the 850 K and 950 K models 
bound the measured luminosity of 2MASS 1207 b, but clearly do not fit the JHK spectral 
shape. Of the A models, the 900 K and 1100 K models bound the measured luminosity of 
2MASS 1207 b, while the 1000 K model gives the best fit. The 1000 K A-model does not 
adequately fit the JHK spectroscopy of 2MASS 1207 b. However, the fit is close enough, 
both in spectral shape, and magnitude so as to suggest that thick clouds might play a roll 
in explaining the under-luminosity of 2MASS 1207 b. The 1000 K model with a modal 
grain size of 30/im fits the K-band spectrum of 2MASS 1207 b, although it overestimates 
the absorption of CO at 2.3 fim and does not fit the H-band spectrum. The 1000 K model 
with a modal grain size of lOO/im fits the H-band spectrum of 2MASS 1207 b, although it 
does not quite match the sharpness of the 2MASS 1207 b H-band spectrum and does not fit 
the K-band spectrum. 

With a dedicated modeling effort, it might be possible to "tweak" the thick-cloud 1000 
K model to fit the spectroscopy of 2MASS 1207 b. As shown in Figure [H frame (c), a 
modal grain-size of 30/im produces a model that is too red in H-K, and a modal grain- 
size of 100/xm produces a model that is too blue. Adjusting the grain-size distribution could 
address this problem. Additionally, the sharpness of 2MASS 1207 b's H-band spectrum is not 
matched by the models. Adjusting surface gravity might address this problem. Finally, the 
models overestimate absorption at 2.3 /xm due to CO. Non-equilibrium chemistry, including 
the developrnent o f more complex organic molecules {C2H2, C2H4, etc.) as described in 



Zahnle et al.l (|2009[ ). might remove some of the carbon from CO and produce a better match 
of this feature. Changing the parameterization of the cloud prescription could also help 
address some of these issues. 



The adequate fit of the iPatience et al.l ( I2OIOI ) spectra by the scaled 1600 K AMES 



DUSTY model makes it tempting to choose this model as bes t describing 2MASS 1 207 b . 



However, we find that the 1000 K thick cloud atmosphere of iMadhusudhan et al.l ( 1201 ll ) 



with some modifications, may be able to fit 2MASS 1207 b's spe ctrum, without an unp hys- 



ical radius scaling. Using the cooling curve parameterizations of iBurrows et al.l (120011 ) and 



2MASS 1207 b's assumed age of 5-12 Myr, a T^g = lOOOK implies that 2MASS 1207 b's 



mass is between bMjup and TMjup- 



6. Conclusions 



We have investigated three hypothetical explanations for the ob served under-luminosit y 
of 2MASS 1207 b: (1) the edge-on disk hypothesis proposed by iMohantv et aD (12003), 
(2) an isotropic dust-shell, motivated by the possible discovery of one around G 196-3 B 
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(IZapatero Osorio et al.ll2010l ). and (3) thick cloud model atmospheres (IMadhusudhan et al. 



20 111 ) that might be able to fit simultaneously explain 2MASS 1207 b's spectral shape and 



luminosity. 

We find the edge-on disk hypothesis unlikely for the following reasons: 

1) Based on modeling and observations, young stars of all masses are expected to exhibit 
variability when occulted by an edge-on disk, as a result of the disk's non-axisymmetric 



structure and hydrostatic fiuctuations. Using data from HST/NICMOS (jSong et al.l 120061 ) 
and VLT/NACO (this work), we see no evidence for strong variability. However, we caution 
that the magnitude and ubiquity of variability in edge-on brown dwarf disks are not currently 
well understood. 

2) The inclination of 2MASS 1207 b has to be very tightly tuned to produce the observed 
under- luminosity effect. At different inclinations, brown dwarfs can be un-extincted, partially 
extincted, or fully extincted by their disks. 2MASS 1207 b falls into the regime of partially 
extincting disks, which our models predict only occur ~6.4% of the time in young brown 
dwarfs. Since 2MASS 1207 b is such a unique system, this low probability is not, in itself, 
enough to completely rule out the edge-on disk hypothesis. 

3) Since the discovery of 2MASS 1207 b and the ensuing discussion of its under- 
luminosity, several more systems have been found that appear to be under-luminous, in- 
cluding HD 203030 B, HN Peg B and HR 8799 bcde. Since these other systems are likely all 
older than 2MASS 1207 b, their gas-rich disks will have dissipated, precluding the geometri- 
cally thick disk geometries necessary to extinct them. Thus at least one other phenomenon 
must be capable of producing the same under- luminosity effect, rendering the edge-on disk 
hypothesis unnecessary. 

While none of our three arguments individually rules out the edge-on disk hypothesis, 
collectively, they strongly suggest that other solutions to 2MASS 1207 b's apparent under- 
luminosity should be considered. One possibility is that an isotropic dust-shell could provide 
the same extinction as an edge-on disk, without the geometric constraints that make an edge- 
on disk unlikely. By enshrouding equilibrium model atmospheres with isotropic spheres of 
dust at different radii, we find the following: 

1) Optically thick dust shells near the surface of the brown dwarf would emit blackbody 
radiation at high enough temperatures so that we would not observe "gray-extinction" in 
the near-infrared. 



2) Optically thick dust shells further from the surface (or at lower temperatures) would 
emit blackbody radiation in the mid-infrared. Our new 8.7/im Gemini/T-ReCS photometry 
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has a 1-0" upper limit of 0.48 mJy and a 3-cr upper limit of 0.92 mJy, which in our simplistic 
models, rules out dust shells at radii less than ~20Rb, or greater than ~250 K. 

Finally, we investigate the possibility that 2MASS 1207 b is not a T^s = 1600 K object, 
despite the fact tha t its spectrum is w ell-fit by a 1600 K AMES-DUSTY model, scaled to 
a l ower luminosity. ICurrie et al.l (l201l[ ) have suggested that thick cloud models (described 
in Madhusudhan et al. 20 111 ) can explain the photometric colors and overall luminosity of 



the HR 8799 planets, which, similar to 2MASS 1207 b, are under-luminous compared to 
AMES- DUSTY models. We a ttemp t to fit 2MASS 1207 b's photometry and spectroscopy 
with the iMadhusudhan et al.l ( 120111 ) atmospheres and find the following: 



1) The 1000 K A-models (thick clouds) are able to reproduce 2MASS 1207 b's low- 
luminosity. The small (amode = 30/im) dust-grain model mostly reproduce 2MASS 1207 
b's K-band spectrum, although it overestimates CO absorption and does not reproduce the 
H-band spectrum. The larger (omode = 100/im) dust-grain model, mostly reproduces 2MASS 
1207 b's H-band spectrum, although it under-estimates the sharpness of the H-band peak 
and does not reproduce the K-band spectrum. 

2) Using Teff = lOOOK and assumin g 2MASS 1207 b's ag e is between 5 and 12 Myr, we 
use the cooling curve scaling relations of iBurrows et al.l (l200l[ ) to estimate that 2MASS 1207 
b's mass is between bMjup and 7Mjup. 

With some "tuning", it seems possible that the thick cloud models will be able to 
explain 2MASS 1207 b's JHK spectrum and photometry, without the need for luminosity 
scaling. Currently, with its trigonometric distance, cluster-membership age and litany of 
spectroscopic and photometric data, 2MASS 1207 b provides the strongest test for the thick 
cloud models. If their fit to the 2MASS 1207 b data is improved, it would signal that 
the under- luminous class of young brown-dwarfs/giant-planets are really cooler, lower- mass 
objects than their spectra imply. It would also demonstrate that young brown-dwarfs/giant- 
planets are capable of suspending dust clouds in their photospheres at cooler temperatures 
than field brown dwarfs. 
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Table 1. Published Photometry/ Spectroscopy of 2MASS 1207 



Reference 


Epoch 


Filter 


A Photometry (mag) 


b Photometry (mag) 


Cutri et al. ('2003') 


May 1999 


J 
H 

Ks 


13.00±0.03 
12.39±0.03 
11.95±0.03 




Javawardhana et al. (2003) 


Apr. 2002 


L' 


11.38±0.10 




Sterzik et al. f2004") 


Jan. 2004 


8.7/xm 
10.4/im 


5.6±1 mjy 
7.5±1 mJy 




Chauvin et al. f2004~) 


Apr. 2004 


H 

Ks 
L' 




18.09±0.21 
16.93±0.11 
15.28±0.14 


Sone et al. f2006~) 


Aug. 2004 
Apr. 2005 


F090M 
FllOM 
F160W 
F090M 
F145M 
F160W 


14.66±0.03 
13.44±0.03 
12.60±0.03 
14.71±0.04 
13.09±0.03 
12.63±0.02 


22.34±0.35 
20.61±0.15 
18.24±0.02 
22.58±0.35 
19.05±0.03 
18.27±0.02 


Mohantv et al. (2007) 


Mar. 2005 
Apr.-Jun. 2005 


J 

HK spectra 


spectra 


20.0±0.2 

spectra 


Riaz et al. (2006) 


Jun. 2005 


IRAC3.6 
IRAC4.5 
IRAC5.8 
IRAC8 
MIPS24 


8.49 mJy 
7.15 mJy 
6.36 mJy 
5.74 mJy 
4.32 mJy 




Morrow et al. (2008) 


Jul. 2005=1 


IRS spectra 


spectra 




Patience et al. (2010) 


Jan.-Feb. 2007 


JHK spectra 


spectra 


spectra 



''IRS spectra of 2MASS 1207 exist for Jul. 2005 and Jul. 2006. iMorrow et"al] ll200d) published the 2005 data, 
but mistakenly listed the date as 2006 (Kevin Luhman, private communication, 2010). 
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Table 2. T-ReCS 8.7 fim Observations of 2MASS 1207 



Date (UT) 


Used On-Source Time (s) 


IQ (%) 


CC (%) 


WV {%) 


BG (%) 


Mar 29, 2008 


304 


70 


50 


50 


50 


Mar 31, 2010 


4445 


70 


50 


50 


Any 


Apr 1, 2010 





70/ Any 


50 


50 


Any 



Note. — Gemini-South Observing Conditions (IQ=Image Quality, CC=Cloud 
Cover, WV=Water Vapor, BG=Background; descriptions available at 
: //www. gemlnl . edu/sciops/telescopes- and- sites/observing- condition- constraints] 



Table 3. Archival NACO Ks photometry of 2MASS 1207 



Date (UT) 


Used On-Source Time (s) 


2MASS 1207 A Ks mag 


2MASS 1207 b Ks mag 


A mag b-A 


2004 Apr. 27 


480 


12.07±0.02 


16.96±0.06 


4.88±0.05 


2005 Feb. 5 


1560 


12.16±0.02 


16.77±0.03 


4.61±0.03 


2005 Mar. 31 


1560 


12.15±0.02 


16.73±0.03 


4.58±0.03 


2009 Jun. 30 


1860 


12.23±0.10 


17.20±0.22 


4.97±0.20 



Note. — Absolute photometry is with respect to the nearby background star, 2MASS 12073400-3932586, which 
has a 2MASS Ks magnitude of 14.86 (and an error of 0.134, which we do not consider above). 



Fig. 1. — Combined Gemini/T-ReCS Si-2 (8.7/im) image of 2MASS 1207, where North is 
up and East is left. 2MASS 1207 A is a t the center of th e two red circles. The near-infrared 
determined position of 2MASS 1207 b (ISong et al.ll2006l ) is at the center of the green circle. 
The halo subtracted region (described in Section 12. ip is the area between the 2 red circles 
(r =0.55" and r =1"), excluding the 2MASS 1207 b aperture, which is marked by the green 
circle (r =0.225"). We find with 50% confidence (i.e. median) that 2MASS 1207 b's flux 
is < 0.26 mJy, 84% confidence (la) that 2MASS 1207 b's flux is < 0.48 mJy and -99.9% 
confidence (3a) that 2MASS 1207 b's flux is < 0.92 mJy 
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l-band Voriability of 2MASS 1207 A 
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Fig. 2.— SMARTS/ANDICAM I-band variability of 2MASS 1207 A. 2MASS 1207 A is 
variable with a modest a ~3%. This means 2MASS 1207 A is a legitimate photometric 
calibrator for 2MASS 1207 b. 
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Fig. 3. — Best-fit fiared and fiat disk models for 2MASS 1207 A. The black, dashed curve is 
the AMES-DUSTY model spectrum for 2MASS 1207 A's photosphere, the red, dot-dashed 
curve is the fiared disk model, the blue, solid curve is the fiat disk model, and the black 



dotted curve is the Spitzer IRS spectrum (IMorrow et al.ll2008l ). Both models adequately fit 
the near- infrared data (< lO/im), but only the fiat model adequately fits the mid-infrared 
data (> lO/xm). Thus, we conclude that the fiat disk model is better for 2M ASS 1207 A, 
and a long with evidence that fiat disk models are better for most brown dwarfs (jSziics et al. 
20101), we use this to justify modeling 2MASS 1207 b's hypothetical disk as fiat. 
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Fig. 4. — Best-fit fiared and fiat disk models for 2MASS 1207 b using four minimum grain 
sizes. The black, dashed curves are the AMES-DUSTY model for 2MASS 1207 b's photo- 
sphere (no extinction), the red, dot-dashed curves are the fiared disk models, and the blue, 
solid curve are the fiat disk models. Note that our 8.7/im upper-limit is drawn with a down- 
ward arrow at the median value connected by a solid line to a horizontal +la upper-limit, 
which is connected by a dashed line to a horizontal +3cr upper-limit. The fiared and fiat 
models are indistinguishable at near-infrared wavelengths (< lO/im), and no data exists at 
longer wavelengths. The smallest minimum grain size (0.1/im) fiat-disk models cannot si- 
multaneously fit the bluest and reddest data. For grain sizes >0.5/im, the models fit the 
data adequately. 
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Fig. 5. — amin = 0.5/im flat disk models for 2MASS 1207 b. For a given model, inclination 
is very tightly constrained to about ~0.3° to fit the observed SED. Small non-axisymmetric 
structures (such as warps, gaps, hydrodynamic fluctutations and dumpiness) in a hypothet- 
ical disk around 2MASS 1207 b, which are analogous to small inclination changes, would 
likely cause variability larger than has been observed (See Section IX2l) . 
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Fig. 6. — The under- luminosity of a brown dwarf extincted by its disk as a function of 
luminosity, shown for our flat disk model, and our flared disk model (the under-luminosity 
of 2MASS 1207 b is marked as a horizontal, dot-dashed line). For low inclinations, the 
disk does not extinct the star. Eventually, the outer edge of the disk moves into the line of 
sight, and extinction rises rapidly. Finally, at very high inclinations, the under-luminosity 
is capped by the scattered light emission of the disk. We term these regions un-extincted, 
partially-extincted, and fully-extincted. That 2MASS 1207 b is in the partially-extincted 
region is, a priori, unlikely, but not statistically refutable. If a disproportionate number of 
objects in 2MASS 1207 b's class are found to be partially-extincted, then there must be 
another source of their under-luminosity. 
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Fig. 7. — Models of 2MASS 1207 b processed through a spherical shell of dust at different 
radii (1-Rb, 5i?b, 10-Rb, 20i?b, 40i?b and 80i?b)- The dashed, black curves are the unattenuated 
Teff = model atmosphere for 2MASS 1207 b and the solid, blue curves are the 

atmosphere viewed through the additional dust shell. Note that our 8.7/im upper-limit is 
drawn with a flat top at +la and a downward arrow at the median value. The models all 
assume 2.5 magnitudes of extinction at J-band, simulating the apparent under-luminosity 
of 2MASS 1207 b. Shells at separations less than 20-Rb are strongly ruled out by our 8.7/im 
+3cr upper-limit. 
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Scaled 1600 K AMES-DUSTY model 
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Fig. 8.— Phot ometry of 2MASS 1207 b compar ed to fa) the b est-fit , scaled AMES- 
DUSTY model jAUard et al boOlL as described in IPatience et aP J201oh . (b)(c)(d) thick 
cloud models (Model-A from iMadhusudhan et al.ll201l[ l at three different temperatures and 
with two different g rain distributions, and (e)( f) intermediate (slightly thinner) cloud mod- 
els (Model-AE from IMadhusudhan et al.ll201ll ) at two different temperatures and with two 
different grain distributions. The 1600 K AMES-DUSTY model is scaled to an unphysical 
object radius of 0.0 5 2 Rp^ (i.e. it is under-luminous or under-sized), whereas the cooler, 
Madhusudhan et al.l ( 1201 ll ) models assume an object size of ~0. 16 Ro), based on the brown- 
dwarf /giant-planet coohng curves of lBurrows et al.l (119971 ). If the IMadhusudhan et al.l (1201 1[ ) 
models are able to fully reproduce the photometry and spectroscopy of 2MASS 1207 b, 
without an unphysical radius scaling, they could explain the whole class of under-luminous 
brown-dwarfs/giant-planets. 
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Fig. 9. — The same models shown in Figure [8] compared to the JHK spectroscopy of 



Patience et ahl ( l2010l ). The 1600 K AMES-DUSTY model reproduces the spectral shape 
of 2MASS 1207 b, but is scaled to a n unphysical object r a dius o f 0.052 Rq (i.e. it is under- 
luminous or under-sized). The cooler iMadhusudhan et al.l (120111 ) models can explain 2MASS 
1207 b's luminosity using a physically motiva ted object radius of ~0.16 R c^, but cannot ex- 



plain its detailed spectral shape. However, the IMadhusudhan et al.l (120111 ) 1000 K A-models, 
shown in (c), are reasonable enough to suggest that modifications to the model might be 
able to explain 2MASS 1207 b's complete appearance. 
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